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This  project  used  parameter  estimation  techniques  to  investigate  dynamic  scaling 
with  an  operational  motorglider  at  the  Test  Pilot  School.  The  Aeromot  AMT  200S  Super 
Ximango  (TG-14A)  was  modified  with  a  data  acquisition  system  to  record  3-axis 
accelerations,  angular  rates,  angle  of  attack  and  sideslip,  Euler  angles,  airspeed,  altitude, 
temperature,  and  control  positions  of  the  rudder,  ailerons,  and  elevators. 

The  primary  objective  was  to  estimate  aircraft  stability  and  control  derivatives  from 
flight  test  data.  The  secondary  objective  was  to  compare  the  derivatives  to  those 
measured  and  dynamically  scaled  during  flight  test  of  a  one-fifth  scale  model.  Flight 
testing  of  the  full  scale  aircraft  was  accomplished  at  Test  Pilot  School  from  14-24  March 
2008.  Flight  test  included  Pitot-static  calibration,  angle  of  attack  and  sideslip  calibration, 
and  pitch  and  yaw-roll  doublets  for  parameter  estimation.  The  stability  and  control 
derivatives  were  estimated  using  a  computational  toolset  called  Systems  IDentification 
Programs  for  AirCraft  (SIDPAC),  a  collection  of  algorithms  for  parameter  estimation. 

The  primary  objective  was  met,  though  some  of  the  data  were  scattered.  The 
secondary  objective  was  not  met.  Some  of  the  parameter  estimates  were  determined  with 
a  high  level  of  confidence,  while  for  others  the  data  were  scattered  and  difficult  to 
interpret.  Nevertheless,  reconstruction  of  time  histories  using  the  final  model  estimates 
produced  data  very  similar  to  the  measurements.  Without  smooth  conditions  normally 
used  for  this  type  of  data  collection  and  programmed  test  inputs  (PTIs),  it  would  be 
difficult  to  improve  the  quality  of  data. 
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INTRODUCTION 

Background 

The  University  of  Arizona  and  Air  Force  Office  of  Scientific  Research  (AFOSR) 
collaborated  in  a  research  program  to  investigate  control  strategies  for  boundary-layer 
separation  of  lifting  surfaces  through  a  combined  approach  of  theory  development, 
numerical  simulations,  laboratory  experiments  and  free-flight  test  experiments. 

Successful  control  of  boundary-layer  separation  had  the  potential  to  lead  to  major 
performance  gains  for  many  Air  Force  applications. 

For  their  research,  the  University  of  Arizona  constructed  a  flight-ready  dynamically 
scaled  1/5  scale  model  of  the  TG-14A  for  wind/water  tunnel  and  free-flight  test 
experiments.  The  University  of  Arizona  wanted  to  compare  experimental  data  from  the 
1/5  scale  to  flight  test  data  of  a  full  size  TG-14A  to  determine  the  “scalability”  of  their 
boundary-layer  flow  control  strategies  over  a  wide  range  of  Reynolds  numbers.  The 
HAVE  MOTO  test  team  used  established  static  stability  and  control  derivative  estimation 
and  flight  test  techniques  to  determine  the  TG-14A’s  stability  and  control  derivatives  for 
comparison  to  experimental  test  results  from  the  1/5  scale  model  as  conducted  by  the 
University  of  Arizona. 

Chronology 

The  USAF  Test  Pilot  School  (TPS),  Class  07B,  conducted  20  flight  tests  totaling  27.1 
hours  at  Edwards  AFB,  California,  from  14  to  25  March  2008. 

Test  Item  Description 


The  TG-14A  (Figure  1)  was  a  two  place  (side  by  side)  T-tail  motorglider  built  by 
Grupo  Aeromot  of  Porto  Alegre,  Brazil,  designed  for  soaring,  pilot  instruction,  and 
training.  The  civilian  designation  for  the  TG-14A  was  the  “AMT  200S  Super  Ximango.” 
In  addition  to  the  standard  aileron,  elevator  and  rudder  control  surfaces,  the  TG-14A  was 
equipped  with  spoilers  on  the  main  wing’s  upper  surface,  retractable  main  landing  gear, 
and  removable  wingtip  extensions  (winglets)  to  increase  the  effective  wingspan.  !  he 
TG-14A  was  powered  by  a  Bombardier-ROTAX  (912S4)  four  cylinder,  liquid  cooled, 
horizontally  opposed  engine.  The  engine  was  rated  at  100  horsepower  at  5800  RPM ,  sea 
level.  Engine  controls  consisted  of  a  throttle,  choke,  and  Propeller  Pitch  Control  Lever 
(PPCL).  For  a  more  detailed  description  of  the  TG-1 4A  refer  to  the  aircraft  flight  manual 
(reference  1). 
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The  TG-14A  used  for  the  HAVE  MOTO  TMP  was  #02-0149  and  was  modified  with 
flight  test  instrumentation  to  provide  time  histories  of  select  data  parameters  required  for 
calculating  stability  and  control  derivatives.  The  instrumentation  consisted  of  a  Cloud 
Cap  Technologies  Piccolo  II  data  collection  system  installed  between  the  aircrew  seats, 
two  additional  air  data  probes  (with  angle  of  attack  and  sideslip  vanes)  installed  at  the 
folding  joints  of  the  two  main  wings  and  string  potentiometer  transducers  installed  on  the 
main  control  surfaces.  The  Piccolo  II  was  linked  to  a  laptop  computer  onboard  the 
aircraft  and  acted  as  the  central  data  acquisition  system  for  all  data  products.  The  Piccolo 
II  also  had  the  capability  to  relay  data  to  a  ground  station  via  a  900MHz  transmitter.  A 
complete  list  of  the  data  parameters  recorded  by  the  Piccolo  II  is  available  in  appendix  A. 


Test  Objectives 

The  overall  objective  of  the  HAVE  MOTO  TMP  was  to  determine  the  stability  and 
control  derivatives  of  the  TG-t4A  for  comparison  to  dynamically  sealed  stability  and 
control  derivatives  from  a  one-fifth  scale  model  of  the  aircraft. 

Objective  1  -  Determine  aircraft  instrumentation  accuracy. 

Objective  2  -  Detennine  stability  and  control  derivatives, 
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Objective  3  -  Compare  stability  and  control  derivatives  to  the  one- fifth  scale  model. 

Objective  4  -  Determine  the  3-dimensional  center  of  gravity  of  the  test  aircraft. 

All  objectives  were  accomplished  with  the  exception  of  objective  3.  No  data  were 
obtained  for  the  scale  model,  and  therefore  no  comparison  was  done. 

Limitations 

The  University  of  Arizona  was  unable  to  fly  the  one-fifth  scale  TG-14  model  as  a 
result  of  Federal  Aviation  Administration  limitations.  Consequently,  objective  3 
(compare  stability  and  control  derivatives  to  the  one-fifth  scale  model)  was  not  fulfilled. 
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TEST  AND  EVALUATION 
Aircraft  Instrumentation  Accuracy 

In  order  to  obtain  accurate  results  for  stability  and  control  derivatives,  the  dynamic 
pressure  and  angle  of  attack  and  sideslip  needed  to  be  accurate.  The  pitot-static  system 
was  calibrated  using  the  cloverlcaf  FTT  and  assuming  zero  total  pressure  error.  Tower 
flyby  and  low  altitude  speed  course  FTTs  were  flown  but  the  data  were  considered  too 
noisy  due  to  poor  conditions.  Angle  of  attack  and  sideslip  vane  corrections  were 
obtained  first  by  correcting  the  vane  angles  to  the  aircraft  body  axis,  then  using  flight 
path  reconstruction  on  data  from  pitch  and  yaw  roll  doublets.  Additionally,  Piccolo  II 
measurements  for  angular  rates  (roll,  pitch,  and  yaw)  and  body  axis  accelerations  were 
checked  for  correct  sense  and  magnitude  using  level  turns,  push  over  pull  up,  and  bank  to 
bank  roll  FTTs  and  comparing  accelerometer  and  gyro  measurements  to  handheld  data. 

A  laser  survey  was  performed  in  hangar  1600  at  Edwards  AFB.  This  survey  used 
precision  laser  measurement  equipment  to  determine  the  relative  angles  of  each  boom 
with  respect  to  the  aircraft  body  axis,  as  well  as  the  position  and  orientation  of  the 
Piccolo  relative  to  the  center  of  gravity.  The  measurements  of  the  Piccolo  can  be  taken 
with  confidence,  however  the  booms  were  extremely  flexible  and  were  disturbed  by  the 
measurement  equipment.  Therefore,  the  numbers  for  boom  offsets  were  not  used. 

Instead,  an  inclinometer  was  used  to  measure  the  angle  of  the  body  axis  and  the  angles  of 
each  boom,  and  the  difference  was  used  for  the  bias.  The  inclinometer  was  less  intrusive 
to  the  measurement  than  the  laser  equipment,  and  no  boom  movement  was  observed 
while  taking  the  measurements.  The  biases  were  subtracted  from  the  measurements  to 
correct  them  to  the  body  axis.  However,  this  was  only  possible  for  angle  of  attack.  Yaw 
offsets  from  the  laser  survey  were  approximately  zero,  and  no  other  method  was  used  to 
estimate  them.  Consequently,  the  sideslip  offsets  were  assumed  to  be  zero.  Therefore, 
they  are  a  potential  source  of  error  in  the  data.  There  was  not  time  to  accomplish  the 
measurement  using  any  other  method.  The  yaw  offsets  for  each  boom  should  be 
determined  using  a  non-intrusive  measurement  technique.  (Rl)1 


Table  1 :  Boom  Angle  Correction  Results 


Inclinometer 

Left  Boom 

Right  Boom 

Angle  of  Attack  Offset 

2.0  deg 

2.9  deg 

Laser  Survey 

Left  Boom 

Right  Boom 

Angle  of  Attack  Offset 

1.5  deg 

2,3  deg 

Angle  of  Sideslip  Offset 

-0.79  deg 

0.29  deg 

Note:  Inclinometer  results  were  used  for  data  reduction,  and  angle  of  sideslip  offset  was 
assumed  zero, 


1  Numerals  preceded  by  an  R  within  parenthesis  at  the  end  of  a  sentence  correspond  to  the  recommendation 
numbers  tabulated  in  the  Conclusions  and  Recommendations  section  of  this  report. 
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Airspeed  and  Altitude  Corrections 

Airspeed  and  altitude  correction  data  were  collected  using  the  Cloverlcaf  flight  test 
technique  (FTT)  as  described  reference  2.  Cioverleaf  test  points  were  flown  over  a  range 
of  airspeeds  from  50  KIAS  to  90  KIAS  and  a  range  of  altitudes  from  5,000  feet  MSL  to 
10,000  feet  MSL.  Position  error  corrections  for  airspeed,  AVpc,  and  altitude,  AHpc,  were 
calculated  using  methods  also  described  in  reference  2.  Data  arc  presented  in  appendix  B 
and  consist  of  plots  of  airspeed  and  altitude  position  error  corrections  as  a  function  of 
instrument  corrected  velocity,  Vic,  for  each  wing  mounted  air  data  probe. 

As  shown  in  each  figure  of  appendix  B  the  position  error  corrections  were  small  in 
magnitude  and  scattered  about  zero.  A  linear  regression  curve  was  fit  to  the  data  in  each 
plot  with  95%  confidence  intervals.  For  the  right  air  data  probe  the  maximum  and 
minimum  AV^  and  AHpc  values  were  +0.78/-0.76  knots  and  +7.0/-6.8  feet  occurring  at 
airspeeds  of  86  and  77  KIAS  respectively.  For  the  left  air  data  probe  the  maximum  and 
minimum  AV^  and  AHpc  values  were  +1.04/-0.42  knots  and  +8. 4/ -4.0  feet  occurring  at 
airspeeds  of  86  and  77  KIAS  respectively. 

In  general,  the  position  error  corrections  were  not  significant  for  calculating  stability 
control  derivatives.  The  right  air  data  probe  tended  to  have  a  slight  negative  bias  while 
the  left  air  data  probed  tended  to  have  slight  positive  bias.  For  calculating  stability 
control  derivatives  the  left  and  right  air  data  were  average  which  further  reduced  the 
effect  of  any  position  errors. 


Accelerations  During  a  Yaw  Roll  Doublet  with  Power  On,  60  Kt AS 


Figure  2:  Measured  Accelerations  vs.  Smoothed  Accelerations  with  Power  On 
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In  order  to  determine  accurate  angle  of  attack  and  sideslip,  the  angles  were  first 
corrected  to  the  body  axis  system  and  then  a  scale  factor  and  bias  for  each  were 
determined  using  flight  path  reconstruction  on  data  from  pitch  and  yaw  roll  doublets. 
This  analysis  was  accomplished  using  an  algorithm  from  SIDPAC  (System  Identification 
Programs  for  Aircraft),  a  MATLAB  ®  program,  to  generate  estimates  for  stability  and 
control  derivatives  (reference  3).  Measured  velocity,  angle  of  attack,  sideslip,  angular 
rates,  Euler  angles,  and  accelerations  were  smoothed  using  an  optimal  global  Fourier 
smoothing  algorithm  from  SIDPAC.  The  smoothed  inputs  were  treated  as  the  actual 
inputs  during  the  rest  of  the  analysis.  An  example  of  smoothed  inputs  is  shown  in  figure 
2.  The  most  noise  in  the  data  occurred  with  the  accelerometers  while  the  engine  was 
running. 

Next,  the  angle  of  attack  and  sideslip  of  each  vane  were  corrected  for  pitch,  roll, 
and  yaw  rates.  Accelerations  were  also  corrected  for  angular  rates  based  on  the  c.g. 
position  relative  to  the  Piccolo.  Finally,  the  algorithm  determined  if  the  accelerations, 
Euler  angles,  and  angular  rates  were  consistent  with  the  air  data  using  kinematic 
relationships.  This  is  also  called  a  data  compatibility  check.  If  the  measurements  were 
perfect,  they  would  be  completely  satisfied  by  the  equations  of  motion.  The  algorithm 
determined  the  scale  factors  and  biases  for  the  measurements  using  an  output  error 
method,  and  reconstructed  the  data  using  these  correction  factors.  The  measurements 
were  corrected  using  the  following  formula: 


corrected  ^bias 


v  measured 


0  +  T 


scalejactor  j 


Table  2:  Angle  of  Attack  Bias  and  Scale  Factors 


Power 

Wingtips 

Airspeed 

(KIAS) 

a  scale  factor 

&  bias 

Mean 

Std  Dev 

Mean 

Std  Dev 

Off 

On 

60 

0.154 

0.077 

0.053 

0.254 

Off 

On 

80 

0,215 

0.016 

0.0240 

0.043 

Off 

Off 

60 

0.206 

0.055 

-0.342 

0.225 

On 

On 

60 

0.133 

0.103 

-0.270 

0.603 

On 

Off 

60 

0.230 

0.018 

-0.003 

0.003 

Note:  Scale  factors  are  unit] 


ess  and  biases  are  in  degrees 


The  flight  path  reconstruction  analysis  was  accomplished  using  doublets  from 
each  of  the  flight  conditions  tested.  For  angle  of  attack,  the  results  from  both  vanes  were 
averaged  to  generate  a  single  angle  of  attack  before  the  data  compatibility  analysis.  The 
angles  of  sideslip  were  kept  separate  and  analyzed  separately  to  determine  individual  bias 
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and  scale  factors  for  each  vane,  and  then  the  two  vanes  were  averaged  after  the  separate 
corrections  were  applied.  Results  are  shown  in  appendix  C.  Tables  2  and  3  contain  the 
means  for  each  data  set  which  were  used  for  subsequent  analysis.  The  bias  for  both  angle 
of  attack  and  sideslip  was  about  zero,  but  the  data  were  noisy.  The  scale  factors  were 
positive  and  easier  to  determine. 


Table  3:  Sideslip  Bias  and  Scale  Factors 


Power 

Wingtips 

Airspeed 

(KIAS) 

Left  Boom 

Right  Boom 

P  scale  factor 

P  bias 

P  scale  factor 

P  bias 

Mean 

Std  Dev 

Mean 

Std  Dev 

Mean 

Std 

Dev 

Mean 

Std  Dev 

Off 

On 

60 

0.202 

0.016 

-0.003 

1.051 

0.127 

0.010 

0.003 

0.529 

Off 

On 

80 

-0.049 

0.060 

-0.302 

1.217 

0.151 

0.007 

-0.084 

0.363 

Off 

Off 

60 

0.243 

0.010 

-1.279 

0.177 

0.116 

0.015 

-0.383 

0.420 

Off 

Off 

80 

0.051 

0.034 

-0.124 

1.371 

0.124 

0.013 

0.515 

0.364 

On 

On 

60 

0.253 

0.026 

-0.255 

0.891 

0.158 

0.012 

-0.126 

0.743 

On 

Off 

60 

0.240 

0.038 

-0.534 

0.590 

0.160 

0.003 

-0.288 

0.366 

Note:  Scale  factors  are  unitless  and  biases  are  in  degrees 


Stability  and  Control  Derivatives 

Three  static  stability  derivatives,  three  control  power  derivatives,  and  three  damping 
derivatives  were  determined  using  parameter  estimation  techniques  on  data  from  pitch 
and  yaw  roll  doublets.  There  were  many  contributors  to  noise  which  induced  scatter  in 
the  data.  The  lateral  data  had  less  scatter  overall  than  the  longitudinal  data,  and  some  of 
the  derivatives  had  very  good  estimates  while  others  were  questionable.  At  the  end  of  the 
analysis,  the  final  model  was  used  to  predict  longitudinal  and  lateral-directional  responses 
and  determine  how  close  the  model  matched  the  flight  test  data.  The  characteristic 
responses  were  the  same,  although  there  were  some  unexplained  differences. 

All  test  points  were  flown  with  the  stick  and  rudders  fixed  during  the  free  response. 
Rudder  was  defined  as  positive  trailing  edge  left.  Aileron  deflection  was  defined  as  one 
half  of  the  left  minus  the  right  aileron  deflection,  where  trailing  edge  down  was  positive 
deflection.  Elevator  was  defined  as  positive  trailing  edge  down.  Two  different  values  of 
dynamic  pressure  were  tested  by  flying  two  different  values  of  calibrated  airspeed,  60 
KIAS  and  80  KIAS.  Flying  constant  calibrated  airspeed  kept  dynamic  pressure  constant 
as  altitude  changed,  Data  were  taken  with  wingtips  on  and  off,  and  also  with  the  engine 
on  and  off.  The  majority  of  the  test  points  were  flown  at  60  KIAS  with  wingtips  on, 
which  was  the  primary  test  condition  where  the  most  confidence  in  the  data  was  desired. 
Maneuvers  with  power  off  were  flown  in  a  steady  glide,  and  power  on  points  were 
obtained  with  maximum  power  while  climbing  back  up  to  the  top  of  the  data  band.  Data 
were  collected  from  5,000  to  10,000  ft  pressure  altitude.  Appendix  A  lists  the  parameters 
recorded  by  the  Piccolo. 
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Measurements  were  smoothed  using  an  algorithm  in  SIDPAC  and  corrected  for  the 
c.g.  offset,  boom  offsets,  and  angular  rates.  To  determine  the  position  of  the  Piccolo  with 
respect  to  the  c.g.,  a  laser  survey  was  accomplished.  The  reference  datum  for  all 
measurements  was  the  tip  of  the  propeller  spinner,  and  the  aircraft  level  reference  was 
determined  from  the  flight  manual.  Results  are  in  table  4. 

After  correcting  for  c.g.  offsets,  bias  and  scale  factor  corrections  were  applied  to  the 
angle  of  attack,  sideslip,  and  acceleration  measurements.  The  stability  and  control 
derivatives  were  estimated  for  each  maneuver  using  an  output  error  algorithm.  A 
confidence  interval  was  generated  on  the  means  of  each  parameter  for  each  condition. 
The  results  for  each  derivative  at  each  condition  are  shown  in  appendix  F.  Standard 
Cramer-Rao  boundaries  are  shown  without  any  correction  factors  applied.  These  values 
are  typical  for  a  high  aspect  ratio,  lightweight  glider,  and  were  successful  in  modeling  the 
actual  outputs  from  the  flight  data. 


Table  4:  Laser  Survey  Results  for  Piccolo  II  Placement 


Piccolo  Position  Offsets 

X 

-23.36  in 

Y 

3.21  in 

Z 

-9.34  in 

Piccolo  Angular  Offsets 

Roll 

1.35  deg 

Pitch 

-0.79  deg 

Yaw 

-0.53  deg 

Note:  these  results  arc  measured  from  the  c.g.  to  the  Piccolo 
Longitudinal  Stability  and  Control  Derivatives 

Longitudinal  data  were  collected  from  pitch  doublets.  The  characteristic  response  of 
the  TG-14  in  pitch  was  highly  damped  with  no  overshoots.  A  typical  maneuver  length 
was  6  seconds  from  pitch  input  to  stable  trimmed  flight.  Several  strategies  were  used  to 
estimate  longitudinal  parameters.  First,  the  algorithm  was  set  up  to  match  only  three 
states  (a,  q,  az)  and  include  only  Cz  and  Cm  parameters  in  the  model.  The  resulting  data 
were  scattered,  and  resulted  in  poor  results  for  V,  0,  and  ax  when  reconstructing  the  data. 
Next,  the  algorithm  was  set  up  to  match  all  six  states  (V,  a,  q,  0,  ax,  az)  using  Cx,  Cz,  and 
Cm  parameters  in  the  model.  For  most  cases,  this  condition  was  too  strict  and  the 
algorithm  could  not  converge  on  a  solution.  The  best  results  were  obtained  by  including 
Cx,  Cm,  and  Cz  parameters  in  the  model  and  matching  only  five  of  the  states  (V,  a,  q,  0, 
az).  Here,  good  results  were  obtained  when  reconstructing  the  data,  with  the  exception  of 
ax,  which  was  generally  small  in  magnitude  and  not  of  concern.  In  some  cases  this 
caused  the  velocity  to  drift  when  reproducing  responses  from  other  inputs.  However,  it 
improved  angle  of  attack,  pitch  rate,  and  g  response  and  the  data  were  more  consistent 
than  previous  attempts.  Therefore,  this  technique  was  applied  to  all  of  the  data  for  the 
final  analysis. 
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Standard,  un corrected  Cramer-Rao  boundaries  were  estimated  for  each  maneuver  and 
are  shown  in  appendix  F.  The  individual  maneuvers  had  tight  tolerances  on  these 
boundaries,  but  between  maneuvers  there  was  some  scatter  on  the  parameter  estimates. 
Appendix  D  contains  sample  time  histories  of  pitch  doublet  responses  and  model 
estimates.  These  estimates  are  based  on  the  parameters  determined  for  that  particular 
maneuver. 

Since  the  flight  operations  took  place  in  March  during  warm  weather  over  the 
lakebeds,  there  was  always  a  small  amount  of  turbulence  present.  Turbulence  introduced 
some  noise  into  the  data,  some  of  which  was  removed  using  smoothing  techniques. 
However,  large  amplitude  inputs  were  used  to  increase  the  signal  to  noise  ratio  and 
maximize  the  effectiveness  of  the  smoothing  techniques.  This  could  have  introduced 
nonlinear  effects  due  to  the  range  of  angle  of  attack  and  sideslip  that  was  experienced 
during  the  maneuvers.  Additionally,  there  was  no  mechanism  for  a  programmed  test 
input  (PTI),  so  every  square  wave  input  was  slightly  different.  Low  weathercock  stability 
and  yaw  damping  made  it  difficult  to  maintain  trimmed  flight  without  sideslip,  so  only 
maneuvers  with  less  than  +/- 1.5  degrees  of  sideslip  were  used.  Pitch  doublets  often 
excited  the  Dutch  roll  mode  in  the  data  because  of  small  aileron  inputs.  These  were  not 
observed  in  flight  but  traces  could  be  seen  in  the  data.  Only  data  with  good  trim 
conditions  and  inputs  were  analyzed,  but  there  were  still  some  variations  within  those 
data  for  sideslip  and  pilot  inputs. 

To  determine  a  final  estimate  for  each  derivative,  the  T  statistic  was  computed  for  a 
95  percent  confidence  interval  on  the  mean  of  each  data  set.  Some  of  the  time  histories 
from  previous  maneuvers  were  reconstructed  using  the  final  results  to  determine  how 
successful  the  parameter  estimation  was.  The  results  are  in  appendix  G.  The  final  model 
has  similar  characteristics  as  the  flight  data,  including  damping  ratio  and  frequency,  but 
the  oscillations  vary  in  amplitude  and  phase. 

Overall,  the  longitudinal  results  were  noisy  and  inconsistent.  The  longitudinal  static 
stability  was  difficult  to  estimate,  and  data  from  engine  on  runs  resulted  in  approximately 
neutral  stability.  Since  the  airplane  did  not  fly  as  if  it  were  neutrally  stable,  these 
estimates  are  incorrect.  Nevertheless,  the  predicted  response  still  follows  that  of  the 
actual  aircraft.  Predicted  responses  were  based  on  a  state  space  model  where  the  A  and  B 
matrices  contained  the  stability  and  control  derivatives.  More  accurate  predictions  could 
be  obtained  with  the  full  set  of  nonlinear  equations  of  motion,  and  may  provide  some 
insight  into  this  situation.  Insight  into  the  results  may  be  gained  by  using  the  full 
nonlinear  set  of  equations  of  motion,  and  by  using  a  nonlinear  aeromodel  for  parameter 
estimation. 

The  pitch  damping  derivatives  also  had  a  lot  of  scatter.  There  does  not  appear  to  be 
significant  difference  between  power  on  and  off,  or  wingtips  on  and  off  conditions,  but  it 
would  be  difficult  to  tell  with  the  scatter  in  the  data.  The  elevator  control  power 
derivative  was  the  most  consistent,  with  the  exception  of  the  power  on,  wingtips  off 
configuration,  where  only  three  runs  were  useable.  It  appeared  to  be  constant  at  all  flight 
conditions.  The  z  force  coefficients,  which  are  related  to  the  lift  curve  of  the  aircraft, 
were  also  determined  and  are  contained  in  appendix  F  with  the  other  data  for 
completeness.  Despite  the  poor  data  quality,  the  reconstructed  outputs  had  the  same 
characteristics  of  the  flight  data,  and  were  successful  at  predicting  the  airplane  response. 
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Lateral-directional  data  were  collected  from  yaw  roll  doublets.  A  two  second  pause 
was  observed  between  the  two  control  inputs.  The  response  was  lightly  damped  and 
exhibited  more  yaw  than  roll.  A  typical  maneuver  was  25  seconds  from  the  initial  rudder 
input  to  stable  trimmed  conditions.  Lateral  directional  derivative  estimates  were 
determined  by  matching  five  states  (p,  r,  p,  q>,  ay)  and  including  Cn,  Q,  and  Cy  parameters 
in  the  model.  The  data  were  more  consistent  than  the  longitudinal  parameters,  both  in  the 
scatter  of  the  data  and  also  in  the  ability  to  reproduce  other  maneuvers  using  different 
initial  conditions  and  inputs.  They  also  varied  less  with  changes  in  conditions. 

Therefore,  all  five  states  were  estimated  and  no  changes  were  made  to  the  initial  output 
error  algorithm. 

Similar  to  the  longitudinal  response,  uncorrectcd  Cramer- Rao  boundaries  were 
estimated  for  each  maneuver  and  are  shown  in  appendix  F.  Appendix  E  contains  sample 
time  histories  of  yaw  roll  doublet  responses  and  model  estimates.  These  estimates  are 
based  on  the  parameters  determined  for  that  particular  maneuver.  In  general,  they  match 
better  than  the  longitudinal  response.  Potential  sources  of  noise  are  the  same  as  in  the 
longitudinal  case.  Additionally,  the  rudder  pedals  were  more  difficult  to  center  than  the 
stick,  making  it  difficult  to  obtain  perfect  trim  conditions  before  and  after  the  maneuver. 
Aileron  inputs  had  a  tendency  to  affect  the  elevator,  and  some  maneuvers  ended  out  of 
trim.  The  tolerance  for  sideslip  was  +/- 1.5  degrees,  and  the  tolerance  for  bank  angle 
(trim)  was  +/-  5  degrees.  Variations  in  these  conditions  as  well  as  input  forms  could 
account  for  scatter  in  the  data. 

To  determine  a  final  estimate  for  each  derivative,  a  95  percent  confidence  interval  on 
the  mean  was  computed  using  the  T  statistic  for  each  data  set.  Some  of  the  previous 
maneuvers  were  reconstructed  using  the  final  model  to  determine  data  quality.  The 
results  are  in  Appendix  G.  Like  the  longitudinal  response,  these  were  developed  using  a 
state  space  model  provided  by  SIDPAC,  and  are  therefore  based  on  linearized  equations 
of  motion.  The  parameters  exhibit  less  variation  between  flight  conditions,  and  match  the 
data  well  in  frequency  and  damping  ratio. 

The  dihedral  (Cip)  had  the  worst  scatter  of  all  of  the  lateral -directional  data.  It  was 
difficult  to  determine  if  there  was  an  effect  between  flight  conditions  because  of  the 
scatter.  The  Cramcr-Rao  boundaries  on  individual  maneuvers  were  worse  than  the  other 
derivatives,  and  the  individual  results  were  scattered.  The  data  suggest  that  the  aircraft 
was  stable,  which  was  expected,  and  the  magnitude  is  similar  to  other  lightweight,  high 
aspect  ratio  gliders.  The  results  for  both  C|P  and  CW  were  consistent  between  maneuvers 
and  also  did  not  change  between  flight  conditions,  which  was  expected.  The  Cramer- 
Rao  boundaries  were  very  tight  on  individual  maneuvers.  The  yaw  damping  derivative 
(Cnr)  was  more  difficult  to  estimate,  and  for  most  conditions  only  three  or  four  useable 
doublets  were  obtained.  Consequently,  there  is  less  confidence  in  these  results.  The 
means  suggest  that  the  parameter  may  be  consistent  over  all  flight  conditions,  but  the 
standard  error  is  too  large  to  make  this  conclusion.  The  rudder  control  power  (Cn5r) 
estimates  were  good,  and  were  mostly  consistent  among  all  flight  conditions.  The  data 
suggest  that  there  may  be  a  power  effect,  although  only  three  points  were  accomplished 
with  the  power  on  and  wingtips  off.  The  weathercock  stability  (Cnp)  was  the  most 
consistent  static  stability  derivative.  The  Cramer-Rao  boundaries  were  small,  and  the 
results  were  consistent.  Sideforce  derivatives  were  estimated  and  the  results  are 
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contained  in  appendix  F  for  completeness.  Overall,  there  were  insufficient  runs  to  make 
significant  conclusions  given  the  scatter  observed  in  the  data. 

3-Dimensional  Center  of  Gravity 

The  test  aircraft  was  positioned  on  three  calibrated  scales  with  each  wheel  resting  on 
a  separate  scale.  Aircrew  were  present  in  the  aircraft  cockpit.  The  scales  supporting  the 
two  main  landing  gear  wheels  were  level  with  each  other  and  the  scale  supporting  the 
rear  wheel  was  positioned  above  or  below  the  horizontal  plane  of  the  main  landing  gear 
as  shown  in  figure  3. 


Figure  3.  TG-14A  CG  Determination 

W  is  total  aircraft  weight. 

H  is  the  vertical  displacement  between  the  surfaces  of  the  two  scales. 

D  is  the  lateral  displacement  between  the  aircraft  main  and  tail  wheels. 

Flm  and  Frm  are  the  measured  forces  from  the  scales  at  the  main  landing  gear  wheels. 

Ft  is  the  measured  force  from  the  scale  at  the  tail  wheel. 

0  is  the  longitudinal  incline  of  the  aircraft. 

The  aircraft  location  of  the  center  of  gravity  in  the  X-Z  plane  was  determined  by  setting 
the  sum  of  the  forces  and  moments  equal  to  zero  at  multiple  inclination  angles.  The 
resulting  system  of  simultaneous  equations  was  solved  to  give  the  location  of  the  center 
of  gravity.  The  center  of  gravity  was  assumed  to  be  at  the  lateral  centerline  of  the 
aircraft,  and  only  x  and  z  positions  were  determined.  The  results  are  shown  in  Table  6. 
The  reference  datum  was  the  tip  of  the  prop  spinner,  and  all  c.g.  positions  were 
determined  with  360  lb  of  pilot  weight.  For  all  aircrew  and  equipment  combinations,  less 
than  20  lb  of  variation  was  observed  in  the  pilot  weight.  The  c.g.  moved  less  than  one 
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inch  in  both  directions  from  full  to  empty  fuel.  These  values  were  used  to  determine  the 
.  relative  position  of  the  Piccolo  with  respect  to  the  center  of  gravity. 


Table  5.  Center  of  Gravity  Location 


Number  of 
Measurements 

Mean 

95%  Confidence  Interval 

X  zero  fuel  (in) 

13 

98.57 

98.15 

98.98 

Z  zero  fuel  (in) 

13 

9.84 

8.12 

11.46 

X  full  fuel  (in) 

12 

97.87 

97.69 

98.05 

Z  full  fuel  (in) 

12 

8.61 

7.75 

9.47 

Note:  Distances  are  measured  with  respect  to  the  tip  of  the  prop  spinner 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  primary  objective  was  met.  Some  of  the  parameter  estimates  were  determined 
with  a  high  level  of  confidence,  while  for  others  the  data  were  scattered  and  difficult  to 
interpret.  Nevertheless,  good  results  were  obtained  using  data  reconstruction  with  the 
mean  values  for  each  parameter  at  each  flight  condition.  Contributing  factors  to  the  noise 
were  lack  of  pure  pilot  inputs,  variations  in  sideslip  during  maneuvers  and  at  trimmed 
conditions,  and  atmospheric  instabilities.  Without  smooth  conditions  normally  used  for 
this  type  of  data  collection  and  programmed  test  inputs  (PTIs),  it  would  be  difficult  to 
improve  the  quality  of  data.  The  c.g.,  air  data  calibration,  and  angle  of  attack  and  sideslip 
calibration  were  accomplished  with  one  exception.  Due  to  lack  of  confidence  in  the  laser 
survey  results  for  the  boom  measurements,  an  inclinometer  was  used  to  recalculate  the 
angle  of  attack  offset,  and  the  sideslip  offsets  were  assumed  zero.  The  yaw  offsets  for 
each  boom  should  be  determined  using  a  non-intrusive  measurement  technique. 

(Rl,  page  5) 

Despite  the  scatter  in  the  data,  reconstruction  of  time  histories  using  the  final  model 
estimates  produced  data  very  similar  to  the  measurements.  The  angle  of  attack  and 
sideslip  matches  were  shifted  in  phase  and  slightly  in  magnitude,  although  they  exhibited 
the  same  frequency  and  damping  characteristics.  However,  the  results  for  static  stability 
did  not  make  sense  with  power  on  for  this  aircraft.  It  is  possible  that  the  linear  state 
space  model  is  not  adequately  addressing  the  effect  of  this  term,  and  a  fully  nonlinear 
model  may  reveal  more  accurate  comparison  of  predicted  to  flight  test  data. 

Additionally,  since  large  pilot  inputs  were  used  to  increase  the  signal  to  noise  ratio  and 
account  for  the  effects  of  turbulence,  the  use  of  a  nonlinear  aero  model  may  be  more 
appropriate  for  estimating  the  parameters. 
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APPENDIX  A  -  LIST  OF  TG-14A  DATA  PARAMETERS 


Parameter 

Description 

Sample 

Rate 

Range 

Resolution 

a* 

Body  axis  acceleration 

50Hz 

-lg  to  lg 

0.001  ft/sec2  ~ 

Body  axis  acceleration 

50Hz 

-lg  to  lg 

0.001  Wsec1 

a. 

Body  axis  acceleration 

50Hz 

-2g  to  6g 

0.001  ft/secz 

P 

Roll  rate 

50Hz 

-100  to +100 
degrees  /second 

0.01  dcg/sec 

q 

Pitch  rate 

50Hz 

*30  to  +30  degrees 
/second 

0.01  deg/sec 

r 

Yaw  rate 

50Hz 

-60  to  +60  degrees 
/second 

0.01  deg/sec 

p, 

Total  Pressure 

50Hz 

1000  lb/ft-  to  2500 
lb/ft2 

0.1  lb/ft2 

ps 

Static  Pressure 

50Hz 

1000  lb/ft"  to  2500 

lb/ft2 

0.1  lb/ft2 

6c 

Elevator  Position 

50Hz 

-20  to  +20  degrees 

0.05  deg 

6r 

Rudder  Position 

50Hz 

-20  to  +20  degrees 

0.05  deg 

6ai 

Left  Aileron  Position 

50Hz 

-20  to  +20  degrees 

0.05  deg 

&ar 

Right  Aileron  Position 

50Hz 

-20  to  +20  degrees 

0.05  deg 

ai 

Left  wing  angle  of  attack 

50Hz 

-5  to  1 5  degrees 

0.05  deg 

Right  wing  angle  of 
attack 

50Hz 

-5  to  15  degrees 

0.05  deg 

Pi 

Left  wing  side  slip 

50Hz 

-10  to  +10  degrees 

0.05  deg 

Pr 

Right  wing  sideslip 

50Hz 

-10  to  +10  degrees 

0.05  deg 

Pdl 

Differential  Pressure  Left 
wing 

50Hz 

0-300  lb/ft2 

0.5  lb/ft2 

Pdr 

Differential  Pressure 

Right  wing 

50Hz 

0-300  lb/ft2 

0.5  lb/ft2 

PS] 

Static  Pressure  Left  wing 

50Hz 

1000  lb/ft2  to  2500 

lb/ft2 

1.5  lb/ft2 

Psr 

Static  Pressure  right  wing 

50Hz 

1000  lb/ftz  to  2500 

lb/ft2 

1.5  lb/ft2 

OAT, 

Outside  Air  Temperature 
left  wing 

50Hz 

-20°  C  to  50°C 

0.1  °C 
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Figure  B-l:  Airspeed  Position  Error  Correction  as  Function  of  Airspeed  (Right 

Boom) 
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Figure  B-2:  Airspeed  Position  Error  Correction  as  Function  of  Airspeed  (Left 

Boom) 
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Figure  B-3:  Altitude  Position  Error  Correction  as  Function  of  Airspeed  (Right 

Boom) 
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Figure  B-4:  Altitude  Position  Error  Correction  as  Function  of  Airspeed  (Left 

Boom) 
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Figure  C-l:  Angle  of  Attack  Bias 
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Figure  02:  Angle  of  Attack  Scale  Factor 
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Figure  03:  Right  Wing  Air  Data  Boom  Sideslip  Bias 
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Figure  C-4:  Right  Wing  Air  Data  Boom  Sideslip  Scale  Factor 
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Figure  C-5:  Left  Wing  Air  Data  Boom  Sideslip  Bias 


C-6 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


5C 

© 

© 

r* 

.= 

o 

(R 

5? 


■+-  r-* 
H  r-B 

S  oT 

axi  ^ 

E  £ 


SS  < 

o2 

djj  © 
a  00 
O  <8 

J>  s 

cz* « 

;as 

£  £ 

.  * 

'*-<  tS 

O  © 

1 
a  w 

©  9 


©  s 


m- 

d 


co 

d 


H - H 

CM 

d 


■s  4 

«  re 
PP  Q 

«  4- 

C3  £ 

Q  H 


s  ^  © 
^  ^ 

O  ©  .2 

^  eS  ? 
^  S  a 
*-«  EL  ° 

P3  i 


O 
V 
^  „ 

£  S 


5  o  £ 
^  V  H 


g  8 
gS 

£  tf 

re  re 
re 
Q  Q 


35  s 
re  © 
C  6 
£  S3 
<  £ 
O  S 

5rt  i 

S3  U 
—  <D 

re  *□ 

«■§ 
£  « 
flu 

a  u 

©  re 

™  cr 

3  1 

rs  u 

H-  ^ 

a  h 

<u 

^  <u 

S  1-1 


^  i  </)  TT- 
ifl  t  -  0X1 


afe 


-i — i- 


-  ^ 

S  9  tfl 
9  ^  < 

£  "So  O 
o  c:  ^ 

Cl  ’© 


-  c 

<5  9  - 

°  a| 

f  D0  Q 
g  C  ^ 
G-  2> 


^  £ 

9  A  3 

&  £2 

§£§ 

°"  5 


it 

G 

>L 

a> 

5 

O 

Cl 


^  <c 


do  O 
c  co 


O  ui 

9  i< 

W.  £1-  yy 

0J  "43  * 

gSs 

“•  * 


c 

9 

i/->  <£ 

00  o 
c  o 


o 

O 


CM 

d 

I 


CO 

d 


a 

i 


jopej  0|BOS  djisopis  woog  yen 


Figure  C-6:  Left  Wing  Air  Data  Boom  Sideslip  Scale  Factor 
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Figure  D-l:  Sample  Time  History  for  a  Pitch  Doublet  (Wingtips  on,  Glide,  60  KIAS) 
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Figure  D-2:  Sample  Results  from  a  Pitch  Doublet  (Wingtips  On,  Glide,  60  KIAS) 
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Figure  D-3:  Sample  Results  from  a  Pitch  Doublet  (Wingtips  On,  Power  On,  60 

KIAS) 
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Figure  D-4:  Sample  Results  from  a  Pitch  Doublet  {Wingtips  On,  Glide,  80  KIAS) 
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Figure  E-l :  Sample  Time  History  for  a  Yaw  Roll  Doublet  (Wingtips  on,  Glide,  60 

KIAS) 


E-l 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Output-Error  Modeling 

Maneuver  length  =  27,88  (sec) 


time  (sec) 


Figure  E-2:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  On,  Glide,  60 

KIAS) 
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Figure  E-3:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  On,  Glide,  60 

KIAS) 
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Figure  E-4:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  On,  Power  On,  60 

KIAS) 


f  (dps)  ay  (9)  phi  (deg) 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Output-Error  Modeling 

Maneuver  length  -  23,62  (sec) 


- Data 

- Model 


Figure  E-5:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  On,  Power  On,  60 

KIAS) 
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Figure  E-6:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  On,  Glide,  80 

KIAS) 
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Figure  E-7:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  On,  Glide,  80 

KJAS) 


E-7 


r  (dps)  p  (dps)  beta  (deg) 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Output-Error  Modeling 

Maneuver  length  =  22.78  (sec) 


Figure  E-8:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  Off,  Glide,  60 

KIAS) 
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Figure  E-9:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  Off,  Glide,  60 

KIAS) 


E-9 


r  (dps)  p  (dps)  beta  (deg) 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Output-Error  Modeling 


Maneuver  length  =  25,1  (sec) 


time  (sec) 


Figure  E-10:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  Off,  Glide,  80 

KIAS) 
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Figure  E-ll:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  Off,  Glide,  80 

K1AS) 
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Figure  £-12:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  Off,  Power  On,  60 

KIAS) 
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Figure  E-13:  Sample  Results  from  a  Yaw  Roll  Doublet  (Wingtips  Off,  Power  On,  60 

KJAS) 
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Figure  F-1.  CMa  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-2.  Cm„  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  F-3.  Cm«  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  F-4.  Cm0  (Power-Off,  Wingtips-On,  80  KIAS) 
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Figure  F-5.  Cm„  (Power-On,  Wingtips-On,  60  KIAS) 
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Figure  F-6.  Cm0  (Power-On,  Wingtips-Off,  60  KIAS) 
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Figure  F-7.  Cm,  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-8.  Cjviq  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  F-9.  Cm,  (Power-Off,  Wingtips-Off,  60  K1AS) 
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Figure  F-10.  Cjyiq  (Power-Off,  Wingtips-On,  80  KIAS) 
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Figure  F-ll.  C\iq  (Power-On,  Wingtips-On,  60  KIAS) 


F-12 


May  2Q0S 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


(UBjpBJ/ J.)  ib|AI0 


Figure  F-12.  Cmq  (Power-On,  Wingtips-Off,  60  KIAS) 
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Figure  F-13.  CMfic  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-14.  Cm8C  (Power-Off,  Win  grips -On,  60  KIAS) 
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Figure  F-15.  CM6e  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  F- 1 6.  Cms®  (Power-Off,  Wingtips-On,  80  K1AS) 
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Figure  F-17.  Cm5C  (Power-On,  Wingtips-On,  60  KIAS) 
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Figure  F-18.  Cmsc  (Power-On,  Wingtips-Off,  60  K1AS) 
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Figure  F-19.  C„p  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-20.  C„p  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  F-21.  C„p  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  F-22,  C0p  (Power-Off,  Wingtips-On,  80  KIAS) 
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Figure  F-23,  Cnp  (Power-Off,  Wingtips-Off,  80  KIAS) 
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Figure  F-24.  C„p  (Power-On,  Wingtips-On,  60  KIAS) 
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Figure  F-25.  C„p  (Power-On,  Wingtips-Off,  60  KIAS) 
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Figure  F-26.  Cnr  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-27.  C„r  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  F-28.  Cnr  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  F-29.  Cnr  (Power-Off,  Wingtips-On,  80  K1AS) 
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Figure  F-30.  Cnr  (Power-Off,  Wingtips-Off,  80  KIAS) 
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Figure  F-3L  Cnr  (Power-On,  Wingtips-On,  60  KIAS) 
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Figure  F-32.  Cnr  (Power-On,  Wingtips-Off,  60  KIA8) 
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Figure  F-33.  CBsr  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-34.  C„*r  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  F  35.  Cn6r  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  F-36.  C„j>r  (Power-Off,  Wingtips-On,  80  K1AS) 
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Figure  F-37.  C„gr  (Power-Off,  Wingtips-Off,  80  K1AS) 


F-38 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Figure  F-38.  CBgr  (Power-On,  Wingtips-On,  60  KIAS) 
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Figure  F-39.  Cn5r  (Power-On,  Wingtips-Off,  60  KIAS) 


F-40 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Figure  F-40.  Cip  from  all  test  conditions  with  95%  confidence  intervals 


F-41 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


Figure  F-41.  Cip  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  F-42.  C,p  (Power-Off,  Wingtips-Off,  60  K1AS) 


F-43 


May  2008 


Edwards  Air  Force  Base 
Air  Force  Flight  Test  Center 


CO 


(D 

Q 


D 

z 

£= 

13 

DC 


fN 


O 

O 


O 

O 


O 

O 


O 

d 


o 

o 


o 

d 


o 

d 


o 

d 


o 

i 


(ueipej/O  dto 


Figure  F-43.  Ctp  (Power-Off,  Wingtips-On,  80  KIAS) 
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Figure  F-44,  Cip  (Power-Off,  Wingtips-Off,  80  KIAS) 
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Figure  F-45.  Qp  (Power-On,  Wingtips-On,  60  KIAS) 
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Figure  F-46.  Qp  (Power-On,  Wingtips-Off,  60  KIAS) 
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Figure  F-47.  CJp  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-48.  Cjp  (Power-Off,  Wingtips-On,  60  K1AS) 
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Figure  F-49.  Clp  (Power-Off,  Wingtips-Off,  60  K1AS) 
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Figure  F-50.  C!p  (Power-Off,  Wingtips-On,  80  K1AS) 
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Figure  F-51.  CIp  (Power-Off,  Wingtips-Off,  80  KIAS) 
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Figure  F-52.  Clp  (Power-On,  Wingtips-On,  60  KIAS) 
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Figure  F-53.  C,p  (Power-On,  Wingtips-Off,  60  KIAS) 
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Figure  F-54.  C|&a  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-55.  Ci8a  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  F-56.  Cj6a  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  F-57.  Cj8a  (Power-Off,  Wingtips-On,  80  KIAS) 
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Figure  F-58.  C,aa  (Power-Off,  Wingtips-Off,  80  KIAS) 
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Figure  F-59.  Cu}a  (Power-On,  Wing(ips-On,  60  KIAS) 
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Figure  F-60.  C]ea  (Power-On,  Wingtips-Off,  60  KIAS) 
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Figure  F-67.  Cyp  (Power-On,  Wingtips-Off,  60  KIAS) 
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Figure  F-68.  Cz«  from  all  test  conditions  with  95%  confidence  intervals 
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Figure  F-70.  C2a  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  F-72.  Cz„  (Power-On,  Wingtips-Off,  80  KIAS) 
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Figure  G-l.  Predicted  Longitudinal  Response  (Power-Off,  Wingtips-On,  60  KIAS) 
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Figure  G-2.  Predicted  Lateral  Response  (Power-Off,  Wingtips-On,  60  KIAS) 
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Predicted  Longitudinal  Response  for  Final  Model,  Power  Off,  Wingtips  Off,  60  KIAS 
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Figure  G-3.  Predicted  Longitudinal  Response  (Power-Off,  Wingtips-Off,  60  KIAS) 


time  (sec) 


Predicted  Lateral  Directional  Response  for  Final  Modet,  Power  OffT  Wingtips  Off ,  60  KIAS 
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Figure  G-4.  Predicted  Lateral  Response  (Power-Off,  Wingtips-Off,  60  KIAS) 
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Figure  G-5.  Predicted  Longitudinal  Response  (Power-Off,  Wingtips-On,  80  KIAS) 
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Predicted  Lateral  Directional  Response  for  Pina!  Model,  Power  Off,  Wfngtlps  On,  80  KtAS 
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Figure  (■-(>.  Predicted  Lateral  Response  (Power-Off,  Wingtips-On,  80  KIAS) 
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Predicted  Lateral  Directional  Response  for  Final  Model,  Power  Offf  Wingtips  Offt  SO  KiAS 
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Predicted  Longitudinal  Response  for  Ftnai  Mode/,  Power  OnT  Wfngtlps  On,  60  KtAS 
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Figure  G-8.  Predicted  Longitudinal  Response  (Power-On,  Wingtips-On,  60  K1AS) 
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Predicted  Laterai  Directional  Response  for  Final  Mode/,  Power  On,  Wingtips  On,  60  KIAS 
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Predicted  Longitudinal  Response  for  Final  Model,  Power  On,  Wingtips  Off,  SO  KiAS 
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Figure  G-10.  Predicted  Longitudinal  Response  (Power-On,  Wingtips-Off,  60  KIAS) 
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Predicted  Lateral  Directional  Response  for  Final  Model,  Power  Onf  Wingtips  Off,  60  KfAS 
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APPENDIX  H  -  FLIGHT  LOG 


Date 

Aircrew 

Test  Points  Flown 

14-Mar 

Marii,  Leigh 

Instrumentation  Checkout,  verification  of  Piccolo  installation.  Tower 

Flyby  at  60  and  90  KIAS;  Low  Altitude  Speed  Course  at  60  and  90  KIAS 

17-Mar 

Mari3T  Leigh 

Cloverleaf  FTT  at  9,000  ft  and  50,  60,  70,  80,  90  KIAS,  Wingtips  On 

17-Mar 

Marii,  Leigh 

Pitch  Doublets  at  10,000-8,000  ft,  60  KIAS,  Wingtips  On,  Power  Off 

17-Mar 

Traynor,  Singh 

Tower  Flyby  at  50,  60,  70,  80,  90  KIAS,  Low  Altitude  Speed  Course  at 

50.  60,  70,  80,  90,  100  KIAS,  Wingtips  On. 

17-Mar 

Traynor,  Singh 

Pitch  Doublets  at  10,000-7,000  ft,  60  KIAS  Power  On  and  Off,  Wingtips 

On 

18-Mar 

Traynor,  Singh 

Yaw  Roll  Doublets  at  10,000-7,000  ft,  60  KIAS  Power  On  and  Off, 

Wingtips  On 

19-Mar 

Traynor,  Singh 

Pitch  and  Yaw  Roll  Doublets,  60  KIAS  Power  On  and  Off,  Wingtips  On 

19-Mar 

Traynor,  Singh 

Pitch  and  Yaw  Roll  Doublets,  60  KIAS  Power  On  and  Off  and  80  KIAS 
Power  Off,  Wingtips  On 

19-Mar 

Marii,  Leigh 

Pitch  and  Yaw  Doublets,  Level,  50,60,70,80,90  KIAS  Power  On,  Wingtips 
On,  9,000  ft. 

19-Mar 

Marii,  Leigh 

Yaw  Roll  Doublets,  60  and  80  KIAS  Power  Off,  60  KIAS  Power  On, 
Wingtips  On. 

21-Mar 

Traynor,  Leigh 

Cloverleaf  FTT,  Low  Altitude  Speed  Course  50,  60,  70,  80,  90,  100  KIAS 
,  9,000  ft,  Wingtips  Off. 

21 -Mar 

Traynor,  Leigh 

Pitch  and  Yaw  Roll  Doublets  Wingtips  Off  60  and  80  KIAS  Power  Off  and 
60  KIAS  Power  On,  Wingtips  Off. 

24-Mar 

Traynor,  Leigh 

Pitch  and  Yaw  Rol!  Doublets  Wingtips  On  60  and  80  KIAS  Power  Off  and 
60  KIAS  Power  On,  Wingtips  On. 

24-Mar 

Traynor,  Marlin 

Cloverleaf  FTT  at  50-90  KIAS  and  5,000-10,000  ft,  Wingtips  On 

24-Mar 

Traynor,  Marlin 

Cloverleaf  FTT  at  50-90  KIAS  and  5,000-10,000  ft,  Wingtips  On 

25-Mar 

Marii,  Marlin 

Pitch  and  Yaw  Roil  Doublets  at  10,000-7,000  ft,  60  KIAS  Power  Off  and 
Wingtips  On 

25-Mar 

Marii,  Marlin 

Pitch  and  Yaw  Roll  Doublets  at  10,000-7,000  ft,  60  KIAS  Power  Off  and 
Wingtips  On 

27-Mar 

Traynor,  Marlin 

Pitch  and  Yaw  Roll  Doublets  at  10,000-7,000  ft,  60  KIAS  Power  Off  and 
Wingtips  On 

28-Mar 

Marii,  Marlin 

Cloverleaf  FTT  at  55-65  KIAS  and  5,000-10,000  ft,  Wingtips  On 

28-Mar 

Marii,  Marlin 

Cloverleaf  FTT  at  55-65  KIAS  and  5,000-10,000  ft,  Wingtips  On 

H-l 
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APPENDIX  I  -  LIST  OF  ABBREVIATIONS,  ACRONYMS,  AND 

SYMBOLS 

Abbreviation  Definition  Units 

AFOSR  Air  Force  Office  of  Scientific  Research 

ax  X  acceleration  ft/s2 

ay  Y  acceleration  ft/s2 

az  Z  acceleration  ft/s2 

Ci  Roll  moment  coefficient 

C|p  Roll  damping  1/rad 

C|(j  Dihedral  1  /rad 

Cisa  Aileron  Control  Power  1/rad 

Cio  Roll  Moment  Bias 

Cm  Pitching  Moment  Coefficient 

Cmq  Pitch  Damping  l /rad 

Cma  Longitudinal  Static  Stability  1/rad 

Cm&e  Elevator  Control  Power  1/rad 

Cmo  Pitching  Moment  Bias 

Cn  Yawing  Moment  Coefficient 

Cnr  Yaw  Damping  1/rad 

Cnj}  Weathercock  Stability  1/rad 

Cngr  Rudder  Control  Power  1/rad 

Cno  Yawing  Moment  Bias 

Cx  Chord  Force  Coefficient 

Cy  Side  Force  Coefficient 

CyP  Change  in  Sideforce  due  to  Roll  Rate  1/rad 

Cyr  Change  in  Sideforce  due  to  Yaw  Rate  1/rad 

Cyp  Change  in  Sideforce  dur  to  Sideslip  1/rad 

Cyo  Side  Force  Bias 

Cy  Normal  Force  Coefficient 

Czq  Change  in  Normal  Force  Due  to  Pitch  Rate  1/rad 

Cya  Normal  Force  due  to  Angle  of  Attack  1/rad 

Czse  Normal  Force  due  to  Elevator  1/rad 

C,.o  Normal  Force  Bias 

FTT  Flight  Test  Technique 

PPCL  Propeller  Pitch  Control  Lever 

PTI  Programmed  Test  Input 

p  Roll  Rate  deg/s 

q  Pitch  Rate  deg/s 

r  Yaw  Rate  deg/s 

SIDPAC  System  Identification  Programs  For  Aircraft 

TMP  Test  Management  Project 

TPS  Test  Pilot  School 

V  Airspeed  ft/s 

a  Angle  of  Attack  deg 

p  Angle  of  Sideslip  deg 
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Abbreviation 

Definition 

Units 

AHpc 

Altitude  Position  Correction 

ft 

AVpc 

Airspeed  Position  Correction 

knots 

0 

Pitch  Attitude 

deg 

<P 

Roll  Angle 

deg 

Heading  Angle 

deg 
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APPENDIX  J  -  DISTRIBUTION  LIST 

For  “onsite”  distribution: 

812  TSS/ENTL 

307  E  Popson  Ave,  Bldg  1400,  Rm  1 10 


Edwards  AFB  CA  93524-6630  3 

USAF  TPS/EDT 
Attn:  Mr  Gary  Aldrich 
220  S  Wolfe  Ave,  Bldg  1220 

Edwards  AFB  CA  93524-6485  2 

USAF  TPS/CS 

Attn:  Ms  Dorothy  Meyer 

220  S  Wolfe  Ave,  Bldg  1220 

Edwards  AFB  CA  93524-6485  2 

Maj  Chad  Rythcr 
PSC  Box  1998 
20  W  Popson 

Edwards  AFB  CA  93524  1 

James  Murray 
NASA  DFRC 
M/S  D-2228 

Edwards,  CA  93523  1 

Maj  Rotem  Maril 
23  Eagle  Ct 

Edwards  AFB  CA  93523  1 

Maj  Devin  Traynor 
40  Blackbird 

Edwards  AFB  CA  93523  1 

Wg  Cdr  Gyandecp  Singh 
30  Hercules  St 

Edwards  AFB  CA  93523  1 

Capt  Elliott  Leigh 
27  Eagle  Ct 

Edwards  AFB  CA  93523  1 
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Capt  Michael  Marlin 
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For  “offsite”  distribution: 

Dr.  Hermann  F.  Fasel 
The  University  of  Arizona 
Aerospace  and  Mechanical  Engineering 
PO  Box  210119 
Tucson  AZ  85721-01 19 

AFRL/AFOSR 
Attn:  Dr.  John  Schmisseur 
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Ste  325  Rm  3112 
Arlington  VA  22203 

Defense  Technical  Information  Center  (DTIC) 
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ATTN:  Willis  Smith  (DTIC-OCA) 
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